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Abstract 

This paper investigates the interaction of plane incident waves with a wave 
farm in the open ocean. The farm consists of a periodic array of large flap- 
type wave energy converters. A linear inviscid potential-flow model, already 
developed by the authors for a single flap in a channel, is considered. Asymp- 



totic analysis of the wave field allows to obtain new expressions of the re- 



flection, transmission and radiation coefficients of the system. It is shown 
that, unlike a line of heaving buoys, an array of flap-type converters is able 
to exploit resonance of the system transverse modes in order to attain high 
capture factor levels. Relations between the hydrodynamic coefficients are 
derived and applied for optimising the power output of the wave farm. 

Keywords: Wave energy, wave-structure interaction, oscillating wave surge 
converters 



1. Introduction 

Research on wave energy converters (WECs) has concentrated tradition- 
ally on systems of small floating bodies, like offshore heaving buoys (see 
[fj-jij]). However, the seminal theories on WECs that originated from this 
first scientific approach to wave energy extraction in the 1970s, do not cap- 
ture exhaustively the dynamics of the last-generation WECs. The latter are 
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usually large-scale devices designed to be deployed in arrays, some of them 
in the near-shore environment. For example, while studying the dynamics of 
an offshore heaving WEC in a channel, Srokosz (1980) [3| showed that reso- 
nance of the channel sloshing modes is detrimental to the efficiency of power 
absorption. Conversely, in a recent analysis of a large flap-type WEC in a 
channel, Renzi & Dias [|§] noted that the trapping of transverse modes near 
the flap increases the efficiency of the converter. Because of the image effect 
of the channel walls, this fact is also expected to occur in an infinite array 
of flap-type converters. The aim of this work is to discover the dynamics 
of a system of last-generation flap-type WECs and to outline its similari- 
ties and differences with respect to the systems of the first generation. As 
a result of this analysis, an optimisation criterion for an array of flap-type 
WECs is devised, depending on the physical and geometrical parameters of 
the system. 

In Section [2] the behaviour of an array of converters in the open ocean 
is investigated by taking as a reference the theoretical framework of Renzi 
and Dias [5[. The expressions of the free-surface elevation for the diffracted 
and radiated wave field in the fluid domain are derived accordingly. Analysis 
of the wave motion in the far field allows to obtain new formulae for the 
reflection, transmission and radiation coefficients. Various relations between 
the hydrodynamic coefficients are then shown in Section HI Some of these 
relations correspond directly to Srokosz's results ji| for floating bodies of 
symmetric shape in a channel. Some others, on the other hand, incorporate 
specific properties of the wave field generated by the flap-type converter, not 
considered before, and point out the peculiarity of such WEC with respect 
to the converters of the first generation. The analytical model is validated 
against known theories in the small-gap and in the point-absorber limit. In 
Section a parametric analysis is undertaken for optimising the performance 
of the system. It is shown that the maximum capture factor is attained at 
complete trapping of the transverse modes of the array. When complete 
trapping is not possible, partial trapping can still increase the performance 
of the system. Finally, in Section a practical application of an array of 
large flap-type WECs is devised. Comparison with available data obtained 



with a finite-element numerical code is very satisfactory (see Appendix C). 
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Figure 1: Geometry of the array (a) and the reference flap (b) in physical variables. 

2. Mathematical model 

2.1. Theoretical background 

Consider an in-line array of identical flap-type wave energy converters, 
each hinged on a bottom foundation of height d in an ocean of constant depth 
h', as shown in figure HJ Primes denote dimensional quantities. Monochro- 
matic incident waves of amplitude A' T , period T" and frequency u' = 2%/T' 
are incoming from the right and set the flaps into motion, which is converted 
into useful energy by means of generators linked to each device. Since the 
practical applications of such a system are usually in the nearshore jsj , where 
wave fronts are almost parallel to the shoreline because of refraction, normal 
incidence is also assumed. Let w' and b' be the width of each flap and the 
spatial period of the array, respectively. Then the gap between two consecu- 
tive flaps is a! — b' — w' (see again figured]). A Cartesian coordinate system 
is set, with the x' direction orthogonal to the flaps, the y' axis along the flap 
lineup and the z' axis rising from the undisturbed water level z' = 0, positive 
upwards; t' denotes time. Due to periodicity, the origin of the system can be 
set arbitrarily on any flap, which is therefore identified as the reference flap. 
The analysis is performed in the framework of a linear inviscid potential- 
flow theory for small-amplitude oscillations. The velocity potential $' must 
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satisfy the Laplace equation 

V 2 &(x',y > ,z',t>) = (1) 

in the fluid domain. The linearised kinematic-dynamic boundary condition 
on the free surface reads 

#ft, + g&S = 0, z' = 0, (2) 

where g is the acceleration due to gravity and subscripts with commas denote 
differentiation with respect to the relevant variables. Absence of normal flux 
at the bottom yields 

=0, z = -ti. (3) 

Because of normal incidence of the incoming wave field, periodicity of the 
problem requires 

f{x', y'+mb', z', = / V, y', z', t'), m = 0, ±1, ±2, . . . , y' G (-b'/2, b'/2), 

(4) 

where /' indicates any physical quantity associated to the problem and m 
each of the flaps; m = denotes the reference flap. Extension to an obliquely- 
incident wave field can be easily made [see for example 0]. However, since 
flap-type WECs are usually designed to operate under normally-incident 
waves [see 0], only normal incidence will be considered here. Because of 
the periodicity condition (J4]), the solution to the complete problem can be 
obtained by investigating the wave interaction with the reference flap centred 
at the origin, with \y'\ < b' /2. Symmetry of the problem requires also 

= 0, y' = ±b'/2, (5) 

which can be regarded as a no-flux boundary condition on two imaginary 
waveguides at y' = ±b'/2 (see again figure [T]). Let 0'{t') be the angle of 
rotation of the flap, positive if anticlockwise; then the kinematic boundary 
condition on the flap yields 

$^ x , = -0' tl (t')(z' + ti -c')H(z' + ti -c'), x' = ±0, \y'\ < w'/2, (6) 

where the thin-body approximation has been applied The Heaviside step 
function in ([6]) assures absence of flux through the bottom foundation. The 
problem defined above is formally equivalent to that solved by Renzi and 
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Dias [5| for a single converter in a channel. Here the main arguments of the 
theory in j^] are retraced and applied to the array configuration. First, the 
system ([I])-© is non-dimensionalised as follows [see H, eqn. (2.1)] 

(x, y, z, h, w, a, c) = (x', y', z', h', w', a', c')/b', t = \fg]b't', 

$ = {^gUA'Y 1 6 = (b'/A')6', (7) 

where the wave amplitude scale A' <C b' because of the hypothesis of small- 
amplitude oscillations. In expression (EJ), a = (1 — w) G (0, 1) defines the 
aperture of the array. Then time is factored out by setting 

$(x,y,z,t) = &{<l>(x,y,z)e- iut }, 9{t) = 3ft {<de~ iu)t } , (8) 

with to = yjy I 'gu' [see H, eqn. (2.11)]. The global spatial potential 

<P = <P R + <P S 

is the sum of the radiation potential 4> R and the scattering potential <j) s . The 
latter is in turn decomposed into 



b 1 + 



where 



i A 

7 (x, y, z) = ^— cosh k(z + h)e~ lkx (9) 

u> cosh kn 

is the incident wave potential and <ft D is the diffraction potential. <p R and 4> D 
must be both outgoing at large \x\. In , Aj = A'j/A' is the non-dimensional 
amplitude of the incident wave and k is the wavenumber, corresponding to 
the real solution of the dispersion relationship u 2 = ktanhkh. Following 
the method described in Appendix B of [5] , application of the Green integral 
theorem yields two hypersingular integral equations, in terms of the jump 
in radiation and scattering potentials across the plate [see eqn. (BIO)]. 
Those equations are solved by expanding the jumps in potential into series 
of Chebyshev polynomials of even order [for details, see eqn.s (B11)-(B18) 
of Careful treatment of the singularity [see 0, eqn. (B19)] ultimately 
allows to write the potentials in the reference domain \y\ < 1/2 in a new 
semi-analytical form [see @, eqn.s (B24) and (B25)]. The radiation potential 
is 

+oo P +oo 

(j) R {x, y, Z) = Yl Yl <t>npm(x,y,z), (10) 
n=0 p=0 m=—oo 
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where 

J.R I \ iUjV ry , N /* /, 2 \l/2 

</>rvm\ x >V> z ) = — K n xZ n (z)a {2p)n J (1-u) 

Hi \K n Jx 2 + (y - \wu - m) 2 ) 

x C/ 2p (u) V V J -du, (11) 

W x 2 + (y — \wu — m) 2 

being the Hankel function of the first kind and first order. In jUJ, 
V = iuiG is the complex angular velocity of the flap, the subscript m identifies 
the contribution of each single flap, while the subscript p indicates the order 
of the Chebyshev expansion, JJi v being the Chebyshev polynomial of the 
second kind and even order 2p, p = 0,1, P E N. The subscript n identifies 
the contribution of each depth mode 



\/2 cosh K n (z + h) 
(h + u~ 2 sinh 2 K n h) 



Z n(z) = 7^ : . „ . sl /2 > n = 0,l,..., (12) 



where kq = k, while K n = ik n denote the complex solutions of the dispersion 
relationship 

u 2 = — k n temk n h, n = l,2, ... (13) 

Finally, the a^p)n are the complex solutions of a system of linear equations en- 
suring that (J) R satisfies the kinematic condition on the flap [see [El, eqn.s (B22) 
and (B23)]. This system is solved numerically with a collocation scheme, 
therefore the solution (fTTj) is partly numerical. In summary, (p npm ffTTj) indi- 
cates the n-th depth mode, p-th order potential of the wave field radiated by 
the m-th flap (m-th array mode), moving at unison with all the other flaps. 
The diffraction potential is given by 



p 



<f ) D (x,y,z) = J2 E ( 14 ) 

p=0 m=— oo 



where 



[ f 1 

- kxZ (z)/3 2p J (l-u 2 ) 



(£.<*.».*> = " 



m) 2 



H[ 1] (kJx 2 + (y- \wu 

U 2p (u) V V 7 du. (15) 

J x 2 + (y — ~wu — m) 2 
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In the latter, the fap are the complex solutions of a system of linear equations, 
which ensures that <p D satisfies the no-flux condition on the flap [see[i|, eqn.s 
(B22) and (B23)]. Again, (p^ m indicates the p-th order potential diffracted 
by the m-th flap, in the presence of all the other flaps. Note that in <p D (j!5p 
only the 0-th depth mode is present, as required by the solvability of the 
whole radiation-diffraction problem [see Appendix B.2 of La] . Computational 
aspects involved in the numerical evaluation of ffTTl) and ffTol) are detailed in 
52.2 of d. 



2.2. Body motion 

The equation of motion of the reference flap in the frequency domain is 
that of a damped harmonic oscillator [see eqn. (2.33) offl, namely 

[-u 2 (I + fi) + C-icu(is + v pto )]e = F, (16) 

depending on the moment of inertia of the flap / = I'/(pb' 5 ), on the flap 
buoyancy torque C = C / (pgb' 4 ) and on the power take-off (PTO) coefficient 
Vpto = v'pto/ipb'^sfgff), where p is the water density. The latter parameters 
are assumed to be all known. In f JT6|) 

nw <r\ I V"^ K n{h — c) sinh n n h + cosh n n c — cosh n n h 1 
2V2 [ n=0 K 2 n (h + u~ 2 smh 2 K n h) ' J 

is the added inertia torque [see eqn. (2.34) of 0], while 

ttw ^ , , uj \k(h — c) sinh kh + cosh kc — cosh kh] , s 

^=^^{«oo}^ m 18 

2v/2 k 2 (h + bj~ 2 sinh kh) 

and 

vkwAi oj [k(h — c) smhkh + coshfcc — cosh kh] . . 

r — T7o (I") 

2V2 k 2 (h + u~ 2 sinh A;/i) 

denote, respectively, the radiation damping [see eqn. (2.35) of H| and the 
complex exciting torque [see eqn. (2.36) of|5|. If the PTO system is designed 
such that 
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which corresponds to the optimum PTO damping [see eqn. (2.40) offl, then 
the average generated power over a period is 



2 



[c-{i + /iV 2 



uj 2 



+ V 2 + V 



(20) 



Now, the generated power (120]) is maximum under resonant amplification of 
the body motion, which occurs when 



By substitution of the latter expression into (I20p . the optimum power avail- 
able for extraction from each flap is therefore 



1 \F 



2 



P opt - g - ^' ( 22 ) 

which matches the well-known result of Srokosz j^] . The performance of each 
element of the array is assessed quantitatively by using two main factors. The 
amplitude factor 

A F = (/l - C ^ an|e| (23) 

is defined as the ratio between the flap horizontal stroke and the amplitude 
of the incident waves, G being the solution of the equation of motion ([16]) . 
Finally, the capture factor is defined as the ratio between the power extracted 
per unit flap width and the power available per unit crest length 

°' = F«fe- (24) 

where 

co / 2kh \ 

c ° = Tk{ 1 + ^2kh) (25) 

is the group velocity of the incident waves. Since P = P(a) and w = 1 — a, 
the capture factor (I24p depends intrinsically on the aperture of the array. A 
strength point of the method of \^ is that knowing the coefficients cton and 
Po is sufficient to obtain immediately all the physical quantities describing 
the performance of the device (eqn.s [1714241 . without need to evaluate the 
potentials f flO|) and f JT4|) . The wave motion at large distance from the array 
will be now analysed. 
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3. The far field 



In this section the behaviour of the wave field is investigated at large 
distance from the array. First, consider the radiation potential 4>n pm given 
by (jTTj) . For n > the Hankel function in (TTTT) can be rewritten as 




(26) 



where K n denotes the modified Bessel function of the second kind and order 
n [see §8.407 of 0|. Since K\(z) oc e~ z as z — > oo in fl26l) and hence in ffTTj) . 

npm 



the argument of 4>nvm for n > decays exponentially in the far field, so that 
at leading order 

0, |x| -> oo,n > 0. (27) 

This happens since the modes n > physically represent the parasite waves 
generated by the motion of the flaps. These remain trapped near the de- 



npm 



vice and do not contribute to the wave motion in the far field [10J. As a 
consequence, 



P +oo 

lR 



J2 E N^oo- (28) 

p=0 m=— oo 

Now substituting (|TT|) into (1281) . using the asymptotic expression (|A~5|) with 
(X, Y) = (x,y — wu/2) and the integral formulae (IB. 21) . (1B.4h and finally 
developing some straightforward algebra, yields 



iV cosh k(z + h) ^ , ,. , 

^— — -> Ate ±llqkx cos(2gvn/), x ±oo. (29) 



o; cosh fc/i 

q=0 



In the latter expression, 7 g = - {2qTx/k) 2 and g is the largest integer for 
which 7^ is real, while 

ITT 

A± = TjWcua 00 Z (0), (30) 
^ = =F J ™ "3>(0) ^ E «(^)o(-l) P (2p + 1) J2p+l(g?rW) . (31) 

p=0 H 
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In (IHTl) . e q is the Jacobi symbol, while Ji P +\ is the Bessel function of first kind 
and order 2p + 1. Note that the radiation potential in the far field (1291) is the 
sum of a progressive long-crested wave (term q = 0) and several progressive 
short-crested waves (terms < q < q), which correspond to the propagating 
sloshing modes of the equivalent channel configuration of \^ . Expression fl29l) 
is similar in form to (2.24) of j^j (accounting for the various differences in the 
nomenclature), which gives the far- field expression of the radiation potential 
for a floating body, symmetric with respect to the x axis, in a channel. In 
(2.24) of |3j, however, the coefficients are left in a general form, while here 
they are determined explicitly for the flap-type converter. The same steps 
can be repeated to find the far-field expression of the diffraction potential 
<p D ([HD. By substituting ([I5J, flll2]) and (EH in flE} and developing 

the algebra, the diffraction potential in the far field becomes 

^(x,y) ~ + h ) j^R^^cos^qny), x -+ ±oo, (32) 

oj cosh fcn. z — ' 

(7=0 

where 

Z7T 

i? = -—wuP o Z Q (0), (33) 

o 

i? 9 = -\w uZ o (0) e q J2 & P (-l) p (2p + 1) J ^ q ™^ , (34) 

p=o 

Equation fl32|) is similar to (2.25) of [3], in which, however, the R q are left 
in a general form. Note that the calculation of the coefficients and R q is 
straightforward once the linear system for «(2 P )o and /?2 P is solved [see eqn. 
B23of§. 

3.1. The free-surface elevation 

The amplitude of the free surface in the far field is an important pa- 
rameter in order to assess the impact of the array on the wave climate of 
the surrounding area. Given the total potential <&(x,y, t), the free-surface 
elevation is 

C(x,y,t) = -$t\*=o = &{n(x,y)e- iut } , 

where 

r](x,y) = ioj^ 1 + (f) R + (f) D ) (35) 
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is the relevant complex spatial component. Substituting flUJ) for (j) 1 and the 
far- field expressions (1251) and (1521) for R and (fi D , respectively, yields 



r](x,y) 
where 



Aje -ikx + ^ q=0 (AjR q + VA+)e^ kx cos(2qny), x -> +00 
Ej=o(^ T « + Wl-)e-^ cos(2g7ry), x -> 



-00 



(36) 



T = 1 — i? , Tg = —R q . (37) 

Overall, the free-surface elevation is the sum of a long-crested wave (term 
q — 0) and several short-crested waves (terms < q < q), namely the 
propagating transverse modes of the array. Physically, in (136|) the terms A^ 
represent the gth-mode radiation coefficients, R q is the gth-mode reflection 
coefficient and finally T q represents the gth-mode transmission coefficient. 
They enjoy all the general properties of the analogous terms introduced by 
Srokosz [3J for bodies of symmetric shape in a channel. In addition, such 
coefficients have some specific properties, peculiar to flap-type bodies, which 
derive from their analytical structure, as shown in detail in §U Figure [2] 
shows the behaviour of Rq and T against the nondimensional wavenumber 
k for a typical configuration where a = 1/2. The plots in figure [2] compare 



favourably with those of Williams & Crull [ll|, fig. 3] and Porter & Evans 
0, fig. 2], who studied the scattering of incident waves by an array of thin 
screens. Note the spiky behaviour of the coefficients, with spikes occurring at 
the resonant wavenumbers k = 2(g+l)7r, q — 0, 1, . . . of the transverse short- 
crested waves, for which the (q + l)th transverse mode turns from trapped to 
propagating. In the following, the coefficients R , T and A^ will be shown 
to enjoy some interesting properties and to be very useful for determining 
some relations between the hydrodynamic coefficients of the system. 



4. Derivation of relations for an array of flap-type WECs 

In this section, relations are derived for an array of flap-type WECs, 
based on the results obtained in the previous section. Some of these relations 
correspond directly to Srokosz's results [i[ for floating bodies of symmetric 
shape, while some others incorporate specific properties of the wave field ( 136]) 
generated by the array of flap-type converters. In this sense, such expressions 
are new and point out the peculiarity of flap-type WECs with respect to the 
converters of the first generation. 
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Figure 2: Magnitude of the fundamental reflection and transmission coefficients, respec- 
tively Rq (|33|) and To (j3T|) . versus the non-dimensional wavenumber k. In this layout the 
flap width equals the gap size, i.e. a = w = 1/2. 

^.i. Extended Bessho-Newman relation 

First consider .4q and Ro- From (130]) and (133]) it is immediate to get, 
respectively, arg(^4o ) = arg(«oo) — tt/2 and arg(_R ) = arg(/3 ) — vr/2. Since 
arg(«oo) = arg(/3 ) [ see Appendix C of 5], then the complex coefficients Aq 
and Rq must have the same argument, say 5, for which 

A£ = \A£\e i5 , R = \R \e t5 } (38) 

for any wavenumber k. The same relation, but only for k < 2tt, can be also 
derived from the Bessho-Newman relation 



i 




where ()* denotes the complex conjugate. Expression (139]) is obtained by 
applyingGreen's integral theorem to 4> s and (0 — (f> R *) and corresponds to 
(3.2) of [3|], with small variations due to the difference in the nomenclature. 
Note that the Bessho-Newman relation (139]) is a general form valid for any 
floating body, symmetric with respect to the x axis, in a channel (or for an 
infinite array of such bodies). Considering k < 2it, i.e. q = 0, and using the 
identities Aq = -A£ (seeEUJ), with A$ = \A$ \e iS , and T = 1 - Rq (see 137]) . 
( 139]) becomes: 2Rq — 1 = e 2lS , which implies ( 138] ). However, while with the 
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general Bessho-Newman relation (I39p it is possible to obtain ( 13 8 p only in the 
domain k < 2tc, usage of the explicit forms (l3"Uj) and (13"3"1) . respectively for 
Aq and i?o, has allowed to extend ( 138)) to any wavenumber. Furthermore, 
by using (j37j) and (]3"g]h (|35j) yields 

cos5 = i^+r 1 — \ A t R i\ > ( 4 °) 

for any A;. Expression (140]) is a particular form of the Bessho-Newman rela- 
tion, valid for a periodic array of flap-type converters under normally-incident 
waves. Note that for k < 2ir, i.e. q = 0, all the transverse modes are trapped 
near the array and (HU1) reduces to 

cos5 = |.Ro|, k < 2n. (41) 

4-2. Relation between F and A^ (array Haskind relation) 

Consider the complex exciting torque ( I19p and the fundamental radiation 
coefficient (|3T)|) . Isolating the term a 00 from (|3TJ|) and substituting it into (TT9"]) . 
yields after some algebra 

F = ±2AjA±C g . (42) 

According to ( 142|) . the long-crested component of the radiated wave field is 
sufficient to obtain the exciting torque acting on each flap, for any value of 
k. Furthermore, since Aj and C g are real numbers, ( 142|) requires 

F = \F\e iS . (43) 

Expression ( 142]) can be transformed into physical variables via ([7]) , thus giving 

F> = 2pgb'A! I At'C' g . 

The latter is similar in form to the well-known two-dimensional Haskind 



relation [lOj except for the factor 6', which represents the array spacing. Fi- 
nally, note that (14"2~P is an extension to intermediate water depth of Srokosz's 
equation (4.3) in [3]. 

4-3. Relation between F and R 

The relation between the exciting torque and the fundamental reflection 
coefficient can be easily obtained by isolating (3q from ( )33l) . substituting it in 
( JT9l) together with (Tl2l) and developing the algebra, so that 



„ , „ t&nhkh ( ', cosh kc — cosh k h s . 

2A,B,— jr - (h-c + ). (44) 
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Figure 3: Ratio Rq/Aq (|45j) versus non-dimensional wavenumber k (solid line) and approx- 
imate expression for large k (dashed line). Parameters of the system are b' — 91.6 m, 
ti = 10.9 m, c' = 1.5 m. 



Hence the exciting torque acting on each flap is related to the amplitude of 
the long-crested component of the reflected wave field, for any value of k. 

4-4- Relation between Ro and At 

By equating (l4"2j) and (jHJ) it is immediate to obtain 



R kC g 



At tanh kh(h-c+ cosh k k ^~^ kh ) ' 



(45) 



valid for any k. Physically, ( 145|) measures the ratio between the reflective 
capacity of the system as an array of screens and the radiative capacity of 
the system as an array of wavemakers, oscillating at unison. In short waves, 



where the flaps are deemed to be operating 12J, it is roughly k ^> 1, so that 
()45l) becomes 

A ^ (46) 

At ~ 2{h-c-k^Y K > 

as shown in figureEJ Since the ratio ( 1461) is 0(k 1 ^ 2 ), the diffractive phenomena 
occurring in the system dominate over the radiative ones in short waves. This 
suggests that the effects of diffraction are not to be neglected if an accurate 
description of the system dynamics is to be pursued. 

4-5. Relation between v and At 

Consider expression ( fl8l) . which defines the radiation damping v for the 
reference plate. Isolating S {«oo} in ( |30|) and substituting it into ( TT8l) yields, 
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after some algebra, 



„. t t&nhkh /, cosh kc — cosh kh\ , . 

v = 2H {AS} ——^- e+ ks . nhkh j, (47) 

for any k. Incidentally, by isolating Aq in the array Haskind relation 
and substituting it into (14"T|) . the latter becomes 

, , F 1 t&nhkh ( , cosh kc — cosh kh 

3FJ <; — — > — — — [h-c + 



Aj J kC g \ k sinh kh 



which corresponds to expression (C3) of Renzi & Dias |5j]. Note that 
f )45|) and (j!7|) allow to obtain the exciting torque and the radiation damping 
- and consequently the optimum generated power ff22l) - directly from the 
fundamental reflection coefficient Rq. This is a peculiar property of the flap- 
type converter and does not hold in general for converters of different shape. 

The above relations have been used to check the numerical calculations 
in this paper. In order to assess the accuracy of computations for a given 
equation of the form l.h.s = r.h.s., the relative error 



|l.h.s. — r.h.s. 
Ir.h.sl 



(48) 



is defined. For a typical system configuration (see |Appendix CJ ), taking 40 
array modes, 5 depth modes and 5 terms in the Chebyshev expansion, is 
sufficient to obtain a relative error of O(10~ 16 ) in calculating the Haskind 
relation (|42p and O(10~ 15 ) in calculating the remaining relations (j44l - fj4"T|) . 
Hence the method of solution based on the Green's theorem of [5J reveals to 
be fast convergent and very efficient. In the following, the influence of the 
array aperture on the performance of the system is assessed, based on the 
relations found in this section. 

5. Performance evaluation 

Consider the optimum capture factor 
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obtained by substituting the optimum power output (I2"2~j) into (153]) . By re- 
placing F with (1441) . z/ with (147}) . and by performing some algebra, (1491) can 
be rewritten as 

Cf = 1^1 (50) 

F 2(l-a)cos<5' V ; 

where 5 is still the argument of i?o- According to (|50|) . the performance of the 
array depends on the reflection coefficient magnitude and argument, which in 
turn are functions of the array aperture. Hence the solution of the scattering 
problem alone is sufficient to assess the performance of the system via (150]) . 
This result confirms that diffraction effects are fundamental in wave-power 
extraction from flap-type WECs. Therefore, the empiric criterion for which 
"to absorb waves means to generate waves" [4j, valid for small floating bodies 
in the absence of diffraction, does not apply here in full. In the following, 
expression ( 1qTJ|) will be validated against known theories in the small-gap and 
point-absorber limits. Then the maximum capture factor will be assessed. 

5.1. Small-gap limit 

In the limit a — > 0, the flaps become joined to each other and the system 
is two-dimensional. In this case it is Rq — > 1, because of complete reflection 
of the incident wave in the diffraction problem. Then it is straightforward to 
show that (150]) becomes 

Cf -> ~, (51) 

i.e. the capture factor coincides with the classical hydrodynamic efficiency 
for 2D devices 13 ■ 



5.2. Point- absorber limit 

Consider now the limit w = w'/b' — > 0, for fixed array spacing b'. In this 
limit, the wavelength A = 2n/k of the incident wave is much larger than the 
flap width, A ^> w, and the interaction between the flaps is weak. Hence 
the results of the present theory can be compared to those of Budal [2| and 
Srokosz Q for an infinite array of point absorbers. For such a system, the 
efficiency is assessed via the absorption length 

L' = b' = C F w>, (52) 
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Figure 4: Interaction factor (|55|) versus non-dimensional wavelength k for the point- 
absorber approximation. Parameters of the system are b' = 91.6 m, w' = 4.58 m, 
h' = 10.9 m, c' = 1.5 m. The ratio w = w'/b' — 0.05 justifies the use of the point- 
absorber approximation. The vertical dash-dotted lines show the resonant wavenumbers 
k = 2(q + l)ir; the grey dashed line shows the linear behaviour s = fc/4 for k < 2ir. 




which is the ratio between the power captured by the single device and that 
incident per unit wave crest length. For an array of converters the optimum 
absorption length is 

L Q pt = I opt s > (53) 

where l' opt is the optimum absorption length for an isolated body and s is 
an interaction factor 0, IsJ] . When to < A, each flap can be considered as a 
three-dimensional axisymmetric body, whose optimum absorption length for 
given wavelength of the incident wave is 

^ = 4 6 ' (54) 

where £ = 1 for heave and £ = 2 for surge (see [H, [l2|| ) . By substituting 
( 154)) into (|53|) . then the latter into ( 152]) and employing (150]) for the optimum 
capture factor, the interaction factor becomes 



k \Rq\ 



(55) 



2£ cos S ' 

where £ = 2, because in the point-absorber approximation the flap moves es- 



sentially in surge 12] . Figure H] shows the plot of s versus the non-dimensional 
wavenumber k = k'b' for a typical configuration in which w = 0.05. When 
k < 2ir (i.e. A > 1) all the transverse modes are trapped and the behaviour of 
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the curve is linear: s ~ k'b' / (2£). This agrees formally with the results shown 
by Budal j^] and Srokosz [j] for an infinite array of heaving point-absorbers, 
where £ = 1. When k > 2tt, i.e. A < 1, incomplete trapping of the trans- 
verse modes strongly modifies the behaviour of the curve. The interaction 
factor globally decreases, but spikes occur near the resonant wavenumbers 
k = 2(q + l)7r. Note that this dynamics is different from that shown in 
for a system of heaving point absorbers. In the latter, s drops to zero when 
k — > 2(q+l)ir from the left, so that the trapping of the transverse modes has 
a detrimental effect on the performance of the system 0, see fig. 2]. Here, 
instead, resonance produces local maxima of s near k = 2(q + 1)tt (see again 
figure H]) and therefore is beneficial in increasing the optimum efficiency, even 
away from complete trapping. This happens because the resonance of the 
transverse modes enhances the horizontal (surge) actions and drops the ver- 
tical (heave) loads [l3j. Therefore surging WECs benefit the most from the 
resonant mechanisms activating in an array configuration. 

5. 3. Maximum capture factor 

Consider again the optimum capture factor, given by ( I5(jp . Expression 
(HO]) shows that |i2o|/ cos ^ has a unit upper limit (HTj) when all the transverse 
modes are trapped near the array. This situation is the most favourable for 
energy extraction and yields the maximum capture factor attainable by an 
array of oscillating wave energy converters of given aperture a. Substitution 
of gU) into §DD yields 

C °F Pt = = k<2rr (56) 

for the maximum capture factor. Incidentally, note that fl56l) correspond to 
Srokosz's 0] maximum efficiency E' max = wCp ax for a symmetric floating 
body in a channel. Since a < 1, the maximum capture factor fl56l) for the 
array configuration is larger than the well-known limit value of 1/2, obtained 
in the small-gap approximation a — > (see expression l5T|) . Hence the mutual 
interaction between the flaps, which is responsible for the trapping of energy 
near the array in the form of short-crested waves, can increase the capture 
factor of the system (if. 

6. Application to wave energy extraction 

In applications, the incident wave period T' and wavelength A' are known, 
together with the flap width w'. The array aperture a needs to be optimised 
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so that the capture factor is maximum. Expression (1561) for the maximum 
capture factor would suggest to increase the aperture as much as possible, so 
that a — Y 1 and consequently Cp ax — > oo. However, f )56l) is to be regarded 
as a theoretical upper limit. This is motivated by two reasons. First and 
most important, as a — > 1, then w — 1 — a — y 0. Now recall that w = w'/b', 
being w' the width of the single flap, which in practical applications is large. 
As a consequence, in order to have w — > 0, it must be b' — > oo. In this 
limit, A = X'/b' — > 0, k — > oo and expression (156]) for the maximum capture 
factor is no longer valid. Physically, by increasing the array aperture a, the 
spatial period b' increases so much, that the array is no longer able to trap 
all the transverse modes, resulting in more energy leakage. Second, recall 
that expression (I5"6"|) for the theoretical maximum capture factor is obtained 
under the assumption (|2T]) . i.e. that the self-oscillation frequency of each 
flap is tuned to the frequency of the incoming waves. However, flap-type 
converters are usually designed to avoid this eventuality. At body resonance, 
the stroke of the flap would exceed by far the amplitude of the incident wave 



12j. It is then clear that this condition is undesirable and not compatible 
with the power take-off mechanism [lij]. Away from body resonance, the 
maximum values of Cp attained are smaller than those predicted by (|56|) . 
This is due to the presence of the inertial terms at the denominator of P 
f l20|) . which in turn reduce C F f l24|) . 

The mathematical model of §2] is now applied to determine an optimisa- 
tion criterion for the array aperture a, which maximises the power output of 
an array of flap-type WECs. The configuration investigated here is that of 
an infinite array of converters similar to Oyster 800™p. Each converter has 
a width w' = 26 m and is placed upon a foundation of height d = 1.5 m from 
the bottom of the ocean; water depth is h' = 10.9 m. Monochromatic inci- 
dent waves of amplitude A'j = lm, period T' = 7 s (wavelength A' = 62 m), 
representative on average of the wave climate off the west coast of Ireland 



15] , are considered. For these parameters, several different layouts, from 
compact (a = 0.30) to sparse (a = 0.95), are analysed to determine the 
optimum array aperture a op t- In each case the power P' = pA' 2 b' 3 ^ 2 g 3 ^ 2 P ex- 
tracted by a single flap and the relevant capture factor Cp fT2^|) are calculated 
with the mathematical model of §2J Table [T] shows the selected values of a, 
the corresponding value of the array spacing b', the generated power P' in 



Oyster is a trademark of Aquamarine Power Limited. 



19 



a 


0.3 


0.40 


0.50 


a opt — 0.58 


0.70 


0.95 


b'(m) 


37 


43 


52 


62 


87 


520 


P' (kW) 


504 


564 


660 


795 


574 


605 


c F 


0.60 


0.68 


0.79 


0.95 


0.69 


0.73 


£jmax 


0.71 


0.83 


1 


1.19 


(1.67) 


(10) 



Table 1: Array spatial period b', power output P' , capture factor Cf and maximum 
theoretical capture factor C™ ax for an infinite array of flap-type converters similar to 
Oyster 800™. Different apertures are considered, from compact (a = 0.3) to sparse 
(a = 0.95). Calculations are made with the mathematical model of $2] 



kW, the capture factor Cp and the theoretical maximum C™ ax . The largest 
power output and capture factor are attained at the optimum configuration 
a = a op t = 0.58, which corresponds to A' = b' (i.e. k = 2n), the trapping 
wavelength of the first transverse mode [5j. For a < a opt , A' > b' (i.e. k < 2n) 
and all the transverse modes are perfectly trapped (see tabled]). However, 
since a is small, the theoretical maximum (156]) sets a relatively small upper 
limit for Cf- By increasing a, C™ ax increases and so does the actual capture 
factor Cp, until it reaches its maximum at a = a opt . For a > a opt the theoret- 
ical limit C^ ax still increases, while the actual capture factor Cp decreases. 
This happens since in these cases A' < b' (i.e. k > 27r) and complete trapping 
of the transverse modes is not possible, so that fl56]) does not hold in prac- 
tice. Energy leakage associated to the propagating transverse waves lowers 
the power absorption of the array well below the theoretical maximum val- 
ues. In conclusion, the optimum aperture that maximises the capture factor 
is the one for which A' = b' i.e. 



w' 



aopt = 1 - T7 ' ( 57 ) 

which can be used preliminary design formula. 

The theory exposed here reveals to be useful for the optimisation of the 
efficiency of an infinite array of flap-type wave energy converters in incident 
monochromatic waves of given period. Of course, in real seas superposition 
of different wave components must be considered. The power output may 
thus vary, depending on the coupling between the spectrum energy period 



and the torque peak period [15|, |16(. Further analysis is therefore necessary 
to obtain more accurate estimates of wave power generation in random seas. 
Finally, due to real sea bottom conditions, converters in array are likely to 
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be deployed in a staggered configuration and in a finite number. Ongoing 
work is investigating the dynamics of a finite array of staggered converters 
and will be disclosed in the near future. 

7. Conclusions 

A periodic array of flap-type WECs has been analysed in this work by 
using the semi-analytical model of Renzi & Dias (if. Asymptotic analysis in 
the far field has allowed to obtain new expressions for the radiation, reflec- 
tion and transmission coefficients. Relations have been determined between 
the Oth-mode coefficients and the hydrodynamic parameters of the system. 
Some of these relations constitute an extension to intermediate water of the 
previous results obtained by Srokosz [jjj for an array of floating bodies in 
deep water, while some others are peculiar to the flap-type converter. The 
efficiency of the system, evaluated via the capture factor, has been shown to 
depend on the reflection coefficient magnitude and argument, which in turn 
are functions of the array aperture. This result shows that diffraction effects 
are fundamental in wave-power extraction from flap-type WECs. Unlike a 
line of heaving buoys [3], an array of flap-type WECs can exploit the res- 
onance of transverse modes to attain high capture factor levels, even when 
complete trapping of the transverse modes does not occur. The maximum 
capture factor is attained in the regime of complete trapping, for which the 
amount of energy available for extraction is the largest. Given the wave 
period and the flap width, the capture factor can be maximised by varying 
the spacing between the flaps, such that complete trapping of the transverse 
modes occurs. These results have been obtained under the assumptions that 
the fluid is inviscid and the flow is irrotational. Viscous effects and turbulent 
dissipations may reduce the values predicted here, especially near trapping 
frequencies (see (lif). 
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search project "High-end computational modelling for wave energy systems". 
Discussions with Prof. D.V. Evans and Dr X. B. Chen have been illuminating. 
Numerical data provision by Dr G. Bellotti and Mr A. Abdolali is gratefully 
acknowledged. 
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Appendix A. Asymptotic analysis of a summation 

Consider the sum 

+~ H[ 1] (kJX 2 + {Y- m) 2 ) 

S(X,Y)= y \ (A.l) 

^ JX 2 + (Y- m 2 

m=— oo V v / 

where k is a real positive number. To determine the asymptotic behaviour 
of 5* for \X\ — > oo, first consider expressions (2.7) and (2.13) of [12], which 
together give 



m=— oo g=— oo V V ^ / 

Differentiating the latter by using the property H^' = —H^ and substitut- 
ing the result in flA.lj) yields, after some elementary manipulations, 

2 + °° 
5(X, y) = — sign(X) ^ e i7 ? fc|^l e 2^y j ( A 2 ) 

g=— oo 

where 

7? = " (2W^F- (A.3) 
Now note that the argument of the square root in 7 g (IA.3P is positive only if 
\q\ < /c/(27r). Then define q as the largest integer for which 

k . . . 

151 < S - (A.4) 

For |g| < 7 9 is real and the relevant terms in S flA.2j) are oscillating 
functions of X. On the other hand, for \q\ > \q\, 7 9 is purely imaginary and 
the relevant terms in S decay exponentially with X. Hence neglecting the 
evanescent terms in fl A. 2[) . transforming the exponential in Y with the Euler 
formula and developing some straightforward algebra yields 

+™ h[ 1] (k^/x 2 + (y — 777, o 

S(X,Y) = /v2 (v , 2 — 

m _oo \A 2 + (y - m) 2 

x ^e g e^ fc|x| cos(2g7ry), |X| +oo, (A.5) 

where e q = 2 — 5 0(? is the Jacobi symbol and 5 nm the Kronecker symbol, 
n, m G N. 
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Appendix B. Evaluation of an integral 

Consider the integral 



[1 - u 2 ) l/2 U 2p {u) cos 



2qir [y 



wu\ 

TJ 



du, 



(B.l) 



where p and q are integers and \y\ < 1/2. For q = 0, application of the 
property (7.343) of [9] for the Chebyshev polynomials U 2p gives immediately 



l p0 



(1 - u 2 ) l ' 2 U 2p (u) du=^-S p0 . 



(B.2) 



Now consider the case q > 0. Expanding the cosine in (IB.lj) . performing the 
substitution u = cos 9 and using the property 



U 2p (cos( 



sin[(2p + 1)9} 
sin 9 



yields after some algebra 



■k/2 



cos(2qiry) / 2 sin 9 sin[(2p + 1)6*] cos(z cos 9) d9, 



(B.3) 



where z = qirw. Substituting the identity 

2 sin 9 sin[(2p + 1)9} = cos{2p9) - cos[(2p + 2)9} 
in ( IB. 31) . using the property 

ctt/2 



cos 2p9 cos (z cos 9) d9 



[see §3.714 of |9j, where J2 P is the Bessel function of first kind and order 2p, 
and going back to the original variables yields finally 



L pq 



(1 - u 2 fl 2 U 2p {u) 



cos 



x (2 P + l) Jw(g7ra,) , g>0, 



qw 



du = cos(2g7ry)(— l) p 

(B.4) 



where the relation J n -\(z) + J„ + i(z) = 2nJ n (z)/z [see §8.471 of i| has also 
been used. 
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Figure C.5: (a) Amplification factor vs period of the incident waves, (b) Capture factor 
vs period of the incident waves. Parameters of the system are b' = 91.6 m, w 1 — 18 m, 
h' = 10.9 m, d = 1.5 m. The solid line shows the results obtained with the analytical 
model of Section [21 dots show the results of the numerical model. 



Appendix C. Comparison with numerical model 

In this section the mathematical model of £0 is further validated against 
available numerical results. The latter have been obtained with a finite- 
element numerical model developed by the University of Roma Tre (Italy), 



as detailed in [18j. In the numerical model, the array layout is replaced by 
the equivalent configuration of a single plate of width w'/2 on the side of a 
channel of width b'/2. This allows to speed up the calculations without loss 
of physical meaning. The numerical model solves the equation of motion (OQ), 
with boundary conditions ([2]) on the free-surface, ([3]) on the bottom, (JSJ) on 
the channel lateral walls and ([6]) on the flap. A radiation condition, including 
a source term for generating the desired incoming waves, is imposed at the 
open generation boundary, which also allows the waves reflected back by the 
device to leave the computational domain freely. At the end of the flume, a 
sink term is imposed to simulate an open boundary, where the transmitted 
waves leave the domain freely. The flume length is 3 times the wave length to 
assure that there is enough space for waves to develop and then leave the do- 
main. In the geometry chosen for comparison, each flap has width w' = 18 m 
and the ocean has depth h' = 10.9 m. The foundation is d = 1.5 m tall and 
the spatial period of the array is b' = 91.6 m, which corresponds to an array 
aperture a ~ 0.8. The flap thickness is null in the semi-analytical model 
and equal to 1.8 m in the numerical model. Comparisons between the semi- 
analytical model of §|2] and the numerical results are shown in figure IC5I In 
the left panel the amplitude factor A F (T2"3"|) is plotted versus the period of 
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the incident waves for the geometry described above. Agreement between 
analytical and numerical data is very satisfactory. Overall, Ap > 0.5 for 
all periods considered, meaning that the flaps effectively convert the wave 
motion into pitching motion. Figure IC.5I (right panel) shows the behaviour 
of the capture factor Cp f )24|) versus the incident wave period. Comparison 
between analytical and numerical data is very good at small periods, while 
the numerical model predicts larger values with longer waves. This is likely 
to be a thickness-induced effect, which becomes important with larger oscil- 
lations of the flap at larger periods. However, even at large T' the results 
predicted by the two models are still in satisfactory general agreement. Note 
also that the capture factor is Cp > 0.6 in the interval T e [5, 8] s, indicating 
that flap-type WECs are most effective in short waves. 
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